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We develop a linear response theory of solvation of ionic and dipolar solutes in anisotropic, axially 
symmetric polar solvents. The theory is applied to solvation in polar nematic liquid crystals. The 

Q formal theory constructs the solvation response function from projections of the solvent dipolar sus- 
ceptibility on rotational invariants. These projections are obtained from Monte Carlo simulations 
■ of a fluid of dipolar spherocylinders which can exist both in the isotropic and nematic phase. Based 

on the properties of the solvent susceptibility from simulations and the formal solution, we have 
obtained a formula for the solvation free energy which incorporates experimentally available proper- 
ties of nematics and the length of correlation between the dipoles in the liquid crystal. The theory 
provides a quantitative framework for analyzing the steady-state and time-resolved optical spectra 
and makes several experimentally testable predictions. The equilibrium free energy of solvation is 
anisotropic in the nematic phase and is given by a quadratic function of cosine of the angle between 
the solute dipole and the solvent nematic director. The sign of solvation anisotropy is determined 
by the sign of dielectric anisotropy of the solvent: solvation anisotropy is negative in solvents with 
positive dielectric anisotropy and vice versa. The solvation free energy is discontinuous at the point 
of isotropic-nematic phase transition. The amplitude of this discontinuity is strongly affected by 
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the size of the solute becoming less pronounced for larger solutes. The discontinuity itself and the 
. ^ \ magnitude of the splitting of the solvation free energy in the nematic phase are mostly affected by 

microscopic dipolar correlations in the nematic solvent. Illustrative calculations are presented for 
the Stokes shift and Stokes shift correlation function of coumarin-153 in 4-n-pentyl-4'-cyanobiphenyl 
| (5CB) and 4,4-n-heptyl-cyanopiphenyl (7CB) solvents as a function of temperature in both the ne- 

matic and isotropic phase. 
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I. INTRODUCTION 



The problem of polar solvation is one of the oldest problems of Physical Chemistry which yet is still a field of active 
theoretical and experimental research. The calculation of solvation free energy is particularly complex, since it is 
affected by a variety of contributions including the short-range cavity formation energy, medium-range dispersion and 
induction forces (non-polar solvation), and long-range electrostatic interactions (polar solvation). These components 
often compensate and complement each other when solvents of different polarity are considered. As a result, calcula- 
tions of the overall free energy of solvation are challenging and often require phenomenological parametrization. Many 
phenomena encountered in chemistry (spectroscopy, redox reactions, etc.) are, however, less affected by the cavitation 
free energy than by non-polar and polar solvation, while the latter often dominates in polar solvents. Therefore, much 
O effort over the last 80 years, following the work of Bornjl Onsager^ and Kirkwood 3 has focused on the understanding 
and modeling of electrostatic, polar solvation. 
. i-h , The original Born-Onsager idea of calculating the electrostatic solvation free energy as the continuum dielectric 
response to charges of the solute has found broad applications, in particularly to solvation of large molecules often 
encountered in bio-medical research4£ For smaller solutes, formal liquid-state theories, most notably integral equa- 
tion theories, have found broad application. These theories are normally formulated either in terms of site-site^ or 
multipolar interaction^ potentials. The proliferation of computer simulation techniques has helped to clarify many 
microscopic features of solvation as well as to test and refine the formal models. 

Most of the effort in the field of solvation thermodynamics has been focused on the understanding of solvation in 
isotropic solvents. Anisotropic solvation, important for chemical reactivity in biological membranes, surfaces, and 
liquid crystalline solvents, has attracted relatively little attention. Also, from the side of experiment, almost nothing 
is known about thermodynamics of solvation in liquid crystals. There is a very limited evidence on the solvatochromic 
shift from spectroscopy^ and a few solvation dynamics studiesSiifiiii have been reported. Computer experiment on 
solvation in liquid crystals virtually does not exist. Continuum models, representing the effect of solvent anisotropy by 
a tensorial dielectric constant, have been proposediiSii^ These approaches provide a very useful continuum limit since, 
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for instance, the Onsager problem of solvation of a spherical dipole— has an exact analytical solution for continuum 
nematicsi^i^ 

Despite the progress in using dielectric continuum models, a few fundamental problems still need to be resolved. 
First, limits of the applicability of the continuum approximation to solvation in polar nematics have not been estab- 
lished. Liquid crystals are mostly made of bulky elongated molecules, and it is a priori unclear if continuum models 
can be applied to solvation of solutes of size often comparable to the size of the solvent molecules. Second, it is not 
clear if dielectric response of a liquid crystal to the solute electric field can in principle be represented by a single 
quantity, the dielectric constant, in particular close to the isotropic- nematic phase transition. 

The approach we propose in this paper is based on the recently obtained microscopic solution for dipole solvation^ 
The model is based on the assumption that the solute-solvent interaction potential is given by the interaction of 
the solute charges with the solvent dipolar polarization. The solvation free energy is then expressed through the 
polarization autocorrelation function of the pure solvent without any particular assumptions regarding the solvent 
structure. The theory is thus applicable to an arbitrary isotropic dielectric. The goal of this paper is to generalize 
this approach to the case of a solvent with axial symmetry. For solvation in isotropic liquids, two projections of 
the polarization autocorrelation function, longitudinal and transverse, are sufficient to describe the dipolar response. 
Lowering the symmetry of the solvent requires a larger set of projections. We derive a formally exact expression for 
the free energy of ionic and dipolar solvation in Sec.|n] [Eq. (|46() ]. 

The full formulation of the theory requires projections of the polarization correlation function on rotational invari- 
ants. These are obtained here from computer simulations of a fluid of dipolar hard spherocylinders. The application 
of the theory to experiment requires, however, a solution based on the input parameters available from experiment. 
This formulation is given in Sec. IIIII [Eq. (|58|l ] in form of a linear combination of solutions obtained in the limit of 
zero wavenumber (continuum) and infinite wavenumber. The relative contribution of each component depends on the 
correlation length of dipolar fluctuations in the liquid (distinct from the correlation length of the order parameter 
fluctuations in Landau-de Gennes theory of liquid crystals^) . One of the principle results of this study is a very slow 
approach of the solvation free energy to its continuum limit, thus invalidating continuum approaches to solvation of 
small and medium-size solutes. We study the dependence of the free energy of solvation on the angle between the 
solute dipole and nematic director as well as on temperature when crossing the point of isotropic-nematic transition in 
Sec. lIVI The solvation free energy is shown to pass through a discontinuity at the transition temperature and becomes 
anisotropic in the nematic phase. Also, the Stokes shift correlation function changes from a single-exponential decay 
in the isotropic phase to bi-exponential decay in the nematic phase. Our results are summarized in Sec. 



II. THEORY 



The linear response approximation (LRA) provides a solution for the solvation free energy (strictly speaking, the 
chemical potential) fi in terms of the response function x(ri,r 2 ) which gives the dipolar polarization in the point of 
space ri produced by the external electric field E (r 2 ) at the point of space r 2 : 

P(ri)= / x(r 1 ,r 2 )-E (r 2 )dr 2 , (1) 

where subscript "0" for the variables refer to the solute. The solvation free energy is then 

fi = -\Jp(r 1 )-E (r 1 )dr 1 . (2) 

The dependence on two separate positions, instead of ri — r 2 of homogeneous solvents, reflects the inhomogeneous 
nature of the solvent response in the presence of the repulsive core of the solute expelling the solvent from its volume. 
In k-scape, Eq. J2J becomes 



/' = ~ / ^^EoCkiJ.JCkx.kaJ.EoC-ka). C!) 



Here, the Fourier transform of the electric field is taken over the solvent volume fl excluding the space occupied by 
the solute 

Eo(ki) = / E (r)e ik r dr. (4) 

The solute space is made by the van der Waals repulsive cores of its atoms. The radii of the solute atoms exposed to 
the solvent are augmented by the shortest distance to the solvent dipole, which, for cylindrically symmetric molecules, 
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is equal to the radius of the cylindrical part of the molecule. Further, the second-rank tensor \ is 

Xa/3(kl,k2) = J-{SP a (^ 1 )SPp{--k2))o, (5) 

kB-L 

where <5P(k) is the Fourier transform of the fluctuation of the solvent dipolar polarization. 

The LRA solution is independent of the electrostatic field of the solute and the subscript "0" in the angular brackets 
denotes the statistical average taken at the presence of a fictitious solute with the repulsive core of the real solute but 
the electrostatic solute-solvent coupling turned offi 1 ^ 1 ? In a hypothetical case of an infinitely small solute, x IS equal 
to the dipolar susceptibility of pure solvent \ s which depends on only one wavevector: 

X(ki,k 2 ) = <5 kl ,k 2 Xs(ki), ( 6 ) 

where subscript "s" denotes the solvent. 

In the general case, x(k!,k 2 ) is affected by the presence of the solute and depends on two k- vectors. The effect 
of the solute on solvent response can generally be separated into two major contribution. The repulsive core of the 
solute distorts the local density of the solvent around it. The spherically-symmetric solute-solvent pair correlation 
function ho s (r) is then different from the solvent-solvent pair correlation function h ss (r). This density disturbance 
affects the dipolar polarization and, consequently, the response function. Another, by far more significant, effect of 
the solute on the solvent response function is related to the expulsion of the dipolar polarization from the solute 
volume. In continuum models, this effect is responsible for the surface charge at the dielectric cavity and, when the 
cavity does not coincide with the equipotential surface, results in a transverse component in the dielectric response. 
The Maxwell's dielectric displacement D(r) then differs from the external electric field E (r). 

The exclusion of the dipolar polarization from the solute volume is accounted for in Chandler's Gaussian 
approximationiSiSS resulting in the following equation for the k-space response function^ 

X(ki,k 2 ) = <5 kl ,k 2 X,(ki) - x"(ki) • o (ki - k 2 )x s (k 2 ). (7) 

Here, 5ki,k 2 ~ (27r) 3 <5(ki — k 2 ) and #o(k) is the Fourier transform of the step function, which equals to unity inside 
the solute and is zero everywhere else. Further, in Eq. J5J, 

X"(k)=x,(k)-[x s (k)-x'(k)]- 1 , (8) 

where 

X'(k) = fT 1 / dridrax.fa - r 2 )e^~^ (9) 
Jn 

and integration in Eq. © is over the volume O occupied by the solvent. 

The substitution of Eq. Q into Eq. © results in the chemical potential of solvation given by the sum of two 
components: 

fl = flh + Mcorr- (10) 

The first term, /ih, corresponds to the homogeneous response approximation (subscript "h") which assumes that 
correlations of dipolar polarization are not modified by the solute and x(ki,k 2 ) can be approximated by dipolar 
susceptibility of the pure solvent according to Eq. The only modification introduced by the solute is the cutoff of 
the electric field inside the solute [Eq. (0}]: 



dk - 

2 J (27r) 



-V* = o / 7^3 E (k) • x.(k) • Eo(-k). (11) 



In isotropic solvents, the tensor x s is diagonal in a coordinate system with one axis taken along k. Its eigenvalues, 
the longitudinal \s an d transverse xT projections, are quite different in the range ak < 2n due to the long-range 
nature of the dipole-dipole interaction potential (a is the diameter of the solvent molecules) i2i In particular, xJ(0) 
grows as the dielectric constant e in strongly polar solvents, while X?(0) tends to a constant value. Because of the 
mutual orthogonality of the longitudinal and transverse projections, the homogeneous solvation free energy fih splits 
into the longitudinal (L) and transverse (T) parts, each given by a 3D integral in k: 

/Xh = fj£ + (12) 
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FIG. 1: Dipolar solute in a nematic solvent. The laboratory system of coordinate is chosen to align the z-axis wit the nematic 
director n. mo denotes the direction of the solute dipole, k is the wave-vector. /3 is the angle between the dipole moment and 
the long axis of the solvent molecule. 

where 

-^ T ^^/(^ X " ,T(fc)l ^' T(k)|2 - (13) 

Once the transverse component of the field E$ (k) is non-zero, which happens when either the solute shape or the 
solute electric field deviates from spherical symmetry, /i^ grows linearly with the dielectric constant resulting in the 
"transverse catastrophe" . Therefore, the homogeneous approximation gives reasonable results only for spherical ions 
when the transverse component of dipolar response is eliminated by the symmetry. It is the second term, generally 
represented by a 6D integral over the two k- vectors, that corrects (subscript "corr") for the unphysical behavior of 
the transverse response. 

The correction component /i corr can be evaluated exactly using analytical properties of the response function when 
the electric field is known in the analytical form. In case of a dipole solute, the final solution is conveniently formulated 
in terms of fRef. flak 

M=(Xtr)~ 1 [^(0)Mh+X^(0)Mh], (14) 

where 

Xtr = ~(xf(0) + 2 X J(Q)). (15) 

Unfortunately, this scheme is hard to implement for liquid crystals. Since a liquid crystal has its own symmetry axis 
(Fig. ^| , x s needs to be diagonalized for each value of k making the problem rather complex. Expansion of the solvent 
response function in spherical harmonics22*2£ appears to be a more straightforward way to the solution. Our solution 
below is given for the case of spherical ion ("i") and spherical dipole ("d"). These two solutes are characterized by 
the following expressions for the field Fourier transform in Eq. 10} 

Ejj (k) = Amqa k (16) 

and 

EW(k) = -4 7r ^m .D k . (17) 

In Eqs. I|16|) and i|17|) . qo an d mo are the solute charge and dipolar moment, and E>k = 3kk — 1 is the dipolar tensor. 
Further, i?o is the solute radius, R\ = Rq + a/2 is the distance of closest solute-solvent separation, and k = k/|k|. 
Here we also use the standard notation for the spherical Bessel functions ji(x). 
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The correction term in Eq. l|l(-)[> is given by a double k-integral 

/W = \ J ^^EoCkx) • x"(ki) • flb(ki - k 2 )x s (k 2 ) • Eo(-k a ). (18) 

In order co convert it to a computationally tractable 3D integral, we first introduce a direct-space field 

Fo(r) = J ^e- k r E (k) • x"(k). (19) 

Analytic properties of the response function x" in complex fc-plane allows one to reduce the integration over k to 
the angular integral over the directions of kii& For the field of a spherical dipole, F = is constant within the 
solute 

F (d) = -^/^m .D k .x"(fc = 0). (20) 

For the spherical ion (see below), 

F W = F = 0. (21) 

In order to calculate Fo in Eq. we need to obtain x"(k = 0) from the k — value of the solvent dipolar 

susceptibility, which we consider next. 



Continuum limit 



We use the laboratory Cartesian system of coordinates with z-axis parallel to the nematic director n (Fig. 
The continuum limit for solvent response function can be obtained from Maxwell's material equations with axially 



symmetric dielectric constant characterized by longitudinal (ei| = e z ) and transverse (ej_ = e :l 

k a kp{e a - l)(e /3 - 1) 



e y ) components 



24 



^Xs,a/3(k = 0) = (e Q - 1) 8 a p 



+ -e±)(k-n) 2 



(22) 



Note that Xs(k = 0) is an even function of k; therefore, according to Eqs. (|16fl and (|19fl . Fo = for the spherical ion 
[Eq. CD]- 

To proceed with the dipolar solute, we first calculate the integral 



"7— X! DKa-fXs,~,(j{k = 0) 
n 7 

The matrix A is diagonal with the elements: 

A — A 
-fT-xx — -^yy 



e±-l 
87r(eii - ej 



(2e||-e ± + 2)-^(e ± + 2)^ 



(23) 



(24) 



A„ = ~ . P 1 > [ex + 2 - (e„ + 2)t/i\ , (25) 



where 



l 



iarctan y/e^Jel - l/ y / e\\/e± - 1 , e\\ > ej_ 
bi(^g)/(>Vi^i7S)!q<^ <26) 

Note that A xx and A zz , as well as xo,xx and Xo,zz below, have no singularities at e|| = e± because terms in the square 
brackets are proportional to — ej_ at e\\ — e± <C 1. 
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In the Appendix, we prove the relation 

Xo = Xs(k = 0) 
Then % is diagonal with the elements 



■*(* = o) = /^*.(* 



0). 



and 



X0,xx — X0,yy 



XO.: 



- 1 



87r(e|| - 



(2e,| 



47r(eii 



- 1 - (en ~ 1)V] ■ 



When Eqs. (1221 and l|27|) are used in the definition of \" m Eq. (JSJ), the final result for the field F^ becomes 



Ro 



where a stands for x,y, z and 



R a = 



e a - n a (e a + 2) 



In Eq. the so-called depolarization factors are given b}^ 

e||(^-l) 



1)' 



(27) 



(28) 



(29) 



(30) 



(31) 



e±-C|| 
1 - 



(32) 



In the isotropic limit, when en = ex, one gets n x = n y = n z = 1/3. 

The field Fq defined by Eq. (|19() is a generalization of the reaction field, introduced by Onsager for a point dipole, 2 
to an arbitrary configuration of solute charges in a solute of arbitrary shape. We have shown here that this field 
reduces to expected limits in the case of spherical ionic and dipolar solutes. In the former case, the reaction potential 
created by the polar liquid within the cavity is constant, and the reaction field is zero. In the latter case, the field 
is constant and our expression in Eqs. (|30|I - H32() coincides with the reaction field in an axially anisotropic dielectric 
previously derived for a spherical dipole by solving the Poisson equation^ 

The zero reaction field in the case of a spherical ion eliminates the correction term in Eq. i|l(J|) . This means that 
the solvent response is longitudinal and the dielectric displacement D is equal to the external field Eo. The free 
energy of solvation is then fully determined by the homogeneous solvation term [Eq. Qllfl]. In case of a dipole, the 
solvent response includes a transverse component, the dielectric displacement is not equal to the external field, and 
the correction term is necessary: 



1 a/3 



m . 



where for a spherical solute 



9 (k) = 4nR\ 



3 h(kR 1 ) 



kR 1 



(33) 



(34) 



Equations (|llf> and i.Y.'A) give the correct continuum limit (subscript 
replacement of X s (k) with its value at k = 0: 



for the ion and 



Q 

2i?i 



~2R\ 



£ c") for the solvation free energy after the 



(35) 



[R X + (R Z - R x ) cos 2 ] 



(36) 
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for the dipole. In Eq. (|36|) . 9o is the angle between the solute dipolar moment and the director (Fig.^l. In the limit 
en — ► ej_, Eqs. I|35[) and (|36|l reduce to their well-know isotropic counterparts, the Born formula 1 



2-Ri V e 
and the Onsager formula 2 

(d) _ m 2 e- 1 

Mo (38) 

Note that the cavity radius is not specified in continuum models. However, empirical experience suggests using the 
van der Waals radius Rq in place of the radius of closest solute-solvent approach R\ appearing in microscopic solvation 
models. 



B. Microscopic theory 

The dependence on the orientation of the wave-vector in Eqs. and i.V-'A) can be integrated out by expanding 
the solvent dipolar susceptibility % s in spherical harmonics [Eq. (|A4|) ]: 

Xs,mn 2 (k) = ^Xs, ni nA k ) Y l*-n^-n^k), (39) 
I 

where LOk denotes the orientation of k. The spherical components of the solute electric field can be obtained for the 
ionic and dipolar solutes^ 

^1(^=4^^^)^^) (40) 

and 

B$ x (k) = -le^ylmo^^^CCl^n^n^m + n[)Y ltn ,{mo)Y* ni+< (u> k ). (41) 

In Eq. H41|) . C^ZiW; Jii, ri2> n) are the Clebsch-Gordan coefficients. 22 Using the product rule and orthogonality of 
spherical harmonics^ we obtain 



(0 



OO 

4g 2 J2 J dkf (kR 1 ) Xs ^ n2l (k)^ir^ Ci ^ 



n 1 n 2 l q 
? 2 



C(112; m, n', ni 4- n')C(112; n2, — n', n2 — n')C(22l; n\ + rc/, 712 — n', ni + n2)moymo.- T 



(42) 



4 v / 6mg 



In Eq. 1021), 



(d) 



5-^ E y rffcjt(fc-Ri)Xs,nin22(fc)(-l) n2 C , (112;ni,n 2 ,ni + n 2 )m , ni (m R)_ ni , 



: \/— ; Xs,nm 2 i( fc )> ( 43 ) 
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and the relations between the spherical and Cartesian components of the vector m = m/m are 

mo,o = mo, z , 

= -(mo,i + im , y )/V2, (44) 
rho-i = (too : x - im a ^ y )/V2. 

Similarly, for the second-rank tensor, one has 

(rh R) = m , z R z , 

(m R)i = -(rho, x + im ,y)R x /V2, (45) 
(mo-R)-i = (m , x - im 0t y)Rx/V2. 

Similar relations exist for the components of the second-rank tensor ^ s (k)^ 

For weakly polar solvents the Cartesian components of \ s (k) form a diagonal matrix in the laboratory system of 
coordinates specified in Fig.^ Moreover, the components related to axes x and y are almost equal to each others. In 
this approximation the solvation free energy reduces to 

oc 

.2 



M w = _!|o J dk f o{kRl)A{kl 



^ = "Sf / dk ^ kR ^ W + C W cos2 ' m 

n > ' 



o 

oo 

2 



/4&r = fpr / dkjKkRx) [R xX s, XX 2(k) - (2J2,x.,*»a(A:) + flxX.,x»2(A;))coB 5, floJ , 
1 o 

= 

r"corr w 3 

where 

2 

^( fc ) = 2Xs,xxo(k) + Xs.zzo(k) + - (Xs,zz2(k) - Xs,xx2{k)) , 

5 

1 4 

B(k) = 3Xs,zzo(k) + 7Xs,xxo{k) + - {?>Xs,zz2{k) - \® X s,xx2{k)) - - {Xs,zzi{k) - Xs,xxi{k)) , (47) 

1 12 
C{k) = Xs,zzo(k) - Xs,xxo(k) + - (5xs,zz2(k) + lQXs,xx2{k)) + y {Xs,zzi{k) - Xs,xxi{k)) . 

In Eqs. Ij46(l and (|47(l . Xs,aai(k) are the coefficients of expansion of the solvent response function in Legendre polyno- 
mials Pi(cos8k)- 

X s , QCt (k) = ^2xs,aai(k)Pi{cos8 k ). (48) 



III. RESULTS 



The theory developed in the previous section requires static dielectric constants and dipolar susceptibility of the 
nematic solvent as input. Here we obtain these data from Monte Carlo (MC) simulation of hard spherocylinders 
with embedded point dipolesi- 2 ^ This fluid transforms from isotropic to nematic phase with decreasing density^ NVT 
MC simulations of N — 800 hard spherocylinders were carried out in our previous studyi2& The dipole moment m is 
parallel to the cylinder axis and the aspect ratio of the length L of the cylindrical part of the molecule to its diameter 
<j is equal to 5. The thermodynamic state of this fluid is fully defined by two parameters: the reduced dipole moment 
(m*) 2 = m 2 / (fceTcr 3 ) and the packing fraction rj = (ir/6)p<r 3 (l + 3L/2a), where p = N/V is the solvent number 
density. Details of the simulation protocol are given in Ref.l26l 

Simulations of nematics with high magnitudes of the dipole moment m* are hindered by the tendency of neighboring 
dipolcs to orient in a locally antifcrroclcctric order. Combined with the elongated shape of the spherocylinders, 
local antiparallcl alignment of dipoles creates bottlenecks in the system phase space, which are hard to explore by 
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standard simulation techniques. 27 In addition, the fluid of dipolar spherocylinders becomes smectic at (m*) 2 > 6.0 
and 77 = 0.47. Because of the relatively high aspect ratio of the solvent molecules, the dipole moment in the range 
< (to*) 2 < 6.0 gives a relatively small overall density of dipoles and, therefore, a low dielectric constant. As a result, 
the capabilities of simple hard-core models are rather limited in exploring high-polarity nematics. The structure of real 
polar nematic liquids, which can demonstrate rather high dielectric constants^ is mostly determined by dispersion 
site-site interactions between elongated molecules. 




FIG. 2: S Z z,i(k) (upper panel) and S xx ,i{k) (lower panel) projections of dipolar structure factors for the fluid of hard sphero- 
cylinders: 1 = (solid line), I = 2 (dashed line) and I = 4 (dotted line). The dash-dotted lines refer to the Pade approximation 
[Eq. H58(l ] with A = OAa. The structure factors are obtained from NVT MC simulations 26 of 800 solvent molecules with the 
packing density r\ = 0.47 and the dipole moment m 2 / (ksTa 3 ) = 1.0. The fluid is in the nematic phase with the nematic order 
parameter S2 = 0.85. 



For weakly polar nematics, the tensor of the dipolar susceptibility \ s is nearly diagonal in the Cartesian coordinate 
system specified in Fig.^ with almost equal xx and yy projections. Therefore, the formally exact formulas in Eq. 142(1 
can be replaced by the approximate relation in Eq. I(46() . The solvation free energy then depends on six one-dimensional 
projection Xs,uai{k) [Eq. 1)48(1] . The corresponding structure factors of dipolar polarization 

S aa ,i(k) = (4TT/3y) Xs , a pi(k) (49) 
have been obtained here from equilibrium MC configurations of the fluid of dipolar spherocylinders (Fig.JSJli^ In Eq. 

y = 47TTO 2 p/(9fc B T) (50) 

is the standard dipole density parameter of dielectric theories^ The reduced dipole moment and the packing fraction 
of the system have values (to*) 2 = 1 and r\ = 0.47 (the isotropic-nematic phase transition occurs at 771N ~ 0.407). In 
this thermodynamic state, the nematic order parameter S2 is equals to 0.85, and the longitudinal (ey, parallel to the 
director) and transverse (ex, perpendicular to the director) dielectric constants are 1.89 and 1.06, respectively. 

The noise in the structure factors obtained from simulations (Fig. [2J arises from the fluctuations of the nematic 
director in the laboratory system of coordinates attached to the simulation box. This setup is necessary to insure that 
the wavevectors used to calculate S aa ,i(k) are eigenvectors of the periodic replicas of the cubic simulation cell. The 
calculations show that the longitudinal structure factors (S zz ,i) are significantly larger than the transverse structure 
factors (S xx ,i), as expected for nematics with longitudinal dipolar moment. In addition, the magnitudes of projections 
decrease rapidly with increasing index I. 

Exact analytical expressions are available for the polarization structure factors in k — *■ and k — > 00 limits. The 
continuum limit relates S zz j(0) and S X x,i(fy to the anisotropic static dielectric constants through Eqs. 1(221) . 1(48(1 . and 
((49(1 . To calculate the k — ► 00 limit, we note that the Cartesian components of the dipolar structure factors are given 
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by the following expression 

S Q/3 (k) = (3/JV) X ' k r • (51) 

where N is a number of solvent particles. All terms in Eq. I|51|) with nj 7^ vanish at k — > 00 resulting in 

S a p(k -> 00) = (3/JV) ^(TO iiQ m ii/3 ) = 3(to q to^). (52) 

This yields 

S xx {k — > 00) = S WW (A: — > 00) = 1 - S* 2 P 2 (cos/3), 

5 zz (fc^oo) = l + 25 2 F 2 (cos/3), (53) 
— > 00) = S , 2 , J (fc — > 00) = S^fc — > 00) = 0, 

where (3 is the angle between the dipole moment and the long molecular axis of the solvent molecule (Fig. ^| and 
P 2 (x) is the second Legendre polynomial. 

From Eqs. (|46p. (|47|l . I|49|l . 152(1. and l|53|) . we find the expression for the solvation free energy when the k — > 00 
limit is used for the solvent susceptibility: 

and 

I'm = (l + ^ 2 P 2 (cos/3)P 2 (cos0 o )^) . (55) 

The continuum [k — 0, Eqs. (|33|) and and short wave-length [fc — > 00, Eqs. ifS*^ and l|5Sjl] limits are two 

asymptotes for the solvation free energy obtained by setting, respectively, the constant Xs(0) and Xs(oo) values for 
the solvent susceptibility in the fc-integrals in Eqs. I|ll|) and (|33|l . The advantage of these limiting expressions is their 
simplicity and direct connection to experimentally available properties of liquid crystalline solvents. These limits can 
be used to derive a practically useful analytical formula for fi. The largest contribution to \i comes from the region 
of k where the squared spherical Bessel functions ji{kR\) in Eq. (|46[) has a maximum. This is the region around 
k = for the ionic solute and k « 2/i?i for the dipolar solute. We need, therefore, a continuous approximation for 
the structure factors that generates a weighted linear combination of k = and k — > 00 limits for the solvation free 
energy. 

The projections of the solvent susceptibility on spherical harmonic are smooth functions of the wavevector for 
weakly polar nematics. The Cartesian components of the dipolar structure factors can, therefore, be reasonably well 
approximated by Pade forms interpolating between the k — and k — * 00 limits 



1 + Afk 2 



?«*,!(*) = ' :^ 2 ■ (56) 



This formula introduces a new theory parameter, the polarization correlation length A;. The correlation length can 
be extracted from structure factors obtained from computer simulations by fitting the slope of S a pj(k) vs k 2 to the 
k — > expansion of Eq. I|56(l (k < k max ~ 2/cr). Values of A; for the fluid of dipolar hard spherocylinders depending 
on packing fraction r\ and dipolar density y are shown in Fig. [21 (closed points). A; was found to be rather weakly 
dependent on ry, even through the isotropic- nematic phase transition (Fig. |3 upper panel). The dependence on y was 
obtained at fixed packing fraction 77 = 0.432 and (m*) 2 changing from 1.0 to 6.0. 

We compare these results to the correlation lengths extracted from slopes of the longitudinal structure factor S L (k) 
vs k 2 and the inverse structure factor l/S T (k) vs k 2 . The latter definition corresponds to the Ornstein-Zernike-Debye 
plot used for the scattering function of liquids close to the critical temperature when the structure factor is a decaying 
function of km^ On the contrary, the longitudinal structure factor is a rising function of k requiring the direct expansion 
of S L (k) in k 2 . The correlation lengths for longitudinal and transverse dipolar fluctuations in isotropic liquids are 
substantially different. This is because these two projections mix together harmonics of the pair distribution function 
with different index I 



s L (k) = i+^(h ii °(k)+2h ii2 (k)y / 

S T (k) = l + l(h^(k)-h^(k)). 



(57) 
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3: The length of dipolar correlations in the fluid of hard spherocylinders vs the solvent packing density (upper panel, 
2 = 1.0) and vs dipolar density y (lower panel, 77 = 0.47). The data are obtained from NVT MC simulations. Closed 



FIG 

(m*; 

diamonds refer to zz, I — projection, closed squares indicate xx, I — projection. Open points in the lower panel indicate the 
transverse and longitudinal correlation length obtained for a fluid of dipolar hard spheres with constant density pa 3 = 0.8 and 
changing dipole moment. The dashed lines connect the simulation points. The dotted vertical tine in the upper panel indicates 
the density of the isotropic-nematic phase transition, rjis = 0.407. 



where h lmn (k) is the Hankel transform [Eq. l|A3jl]. 

In contrast to longitudinal and transverse projections of isotropic fluids, the projections S a /3,i(k) correspond to 
the same index I. The anisotropy of the nematic phase relative to the director is then taken out to the Legendre 
polynomial Pi(cos9k) [Eq. (|48|l ] . In addition, we found that A; calculated from different harmonics with the same I 
are approximately equal to each other, at least for weak polar solvents. This is why A; in Eq. I|50l) does not include 
Cartesian projections a, (3. The results for / — are presented in Fig. while data for I > do not converge well 
because of large statistical errors. 

2 1 ■ , ■ , ■ , ■ , ■ 1 




2 4 6 8 10 

R/A 



FIG. 4: Non-locality functions [Eq. 1591 1 and reduced solvation free energy fi//j, c vs the solute size. The dashed and solid lines 
refer to the spherical ion and spherical dipole solutes, respectively. 

With the structure factors given by Eq. I|56() . the solvation free energy is a linear combination of continuum and 
large wavevector limits: 

H = + /(iii/A) (/loo - Ho) ■ (58) 

The function /(i?i/A) represents the contribution of the non-local solvent response, influenced by the finite length 
of dipolar correlations, to the solvation thermodynamics. For the cases of spherical ionic (i) and dipolar (d) solutes, 
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this function is given by the following expressions 

/«(z)=0.5 [l-e- 2x ] /x, 



/ w (i) = 1.5 x z ~l + (x + l) e 



2 -2x 



/X 3 



(59) 



(NO 
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FIG. 5: Solvation free energy of ionic (upper panel, in q^/Ri units) and dipolar (lower panel, in m%/R\ units) solutes vs Ri. 
Shown are the microscopic calculation ["m", Eqs. 146II - H481 1. the Pade form ["Pade", Eq. 1581 1. and the continuum limit ["c", 
Eqs. lE3 and gg) ]. 



Non- locality functions /W(i?i/A) and /^(Ri/A) decay monotonically from one to zero with increasing the solute 
size R± (Fig. 0J . The ionic non-locality function decays faster than the dipolar function indicating that ionic solvation 
is better described by continuum approximation than is dipolar solvation. Notice that formula (|58|l applies to solvation 
in strongly polar nematics, because it has been derived without assuming weak polarity of the solvent (the cross-terms 
omitted in Eq. (|46|l are included in fi c ). Figure demonstrates a satisfactory agreement between the approximate 
solution given by Eq. I|58|l and the exact microscopic theory [Eqs. I|46|) - (|48|l ]. The increase of the solvent dipole 
moment leads to an oscillatory character of the response function, but even in this case, the Pade approximation [Eq. 
(|56|l ] gives reasonably accurate estimates for the solvation free energy in isotropic polar solvents^ 

Equation 1581) . combining the limiting continuum and high wavevector values for ionic [Eqs. (|35l) and i|54|l ] and 
dipolar [Eqs. H36[) and l|55l) ] solutes, is the central results of this study. The solvation free energy requires the following 
solvent parameters: dielectric constants eii and e±, the order parameter S2, the polarization correlation length A. In 
addition, the ionic charge or dipole moment, along with the radius R\ = Rq + cr/2, should be supplied. The dielectric 
constants and the order parameter come from experiment. The correlation length A is not experimentally available 
and, for the sake of interpreting the experiment, is, at the moment, a theory parameter requiring fitting to some 
experimental observable. Simulations of the model fluid of hard spherocylinders suggest magnitudes of A of the order 
A ~ 0.3 - 0.4cr. 



IV. THEORY PREDICTIONS 



The present theory allows us to make some specific predictions regarding equilibrium solvation and solvation dy- 
namics. The electrostatic component of solvation free energy can be measured from the steady-state Stokes shift of 
optical lines, whereas solvation dynamics is probed by the Stokes shift correlation function. Our discussion below will 
therefore target these two properties. 



13 



A. Equilibrium solvation 



Our present development highlights several issues important for the understanding of equilibrium solvation in 
axially-symmetric solvents: (1) Solvation anisotropy, i.e. the dependence of the free energy of solvation on the orien- 
tation of the solute dipole relative to nematic director. (2) The effect of crossing the phase transition temperature on 
the solvation thermodynamics. (3) The effect of dipolar correlations on solvation and the applicability of continuum 
models of solvation. 

From Eqs. (|3fc>|> . (|55|l . and (|58p. the solvation free energy of a dipole, ^jS d \ is a quadratic function of cos£?o, where 
6q is the angle between the solute dipole and nematic director (Fig. QJ. Figure shows that the solvation anisotropy 



(60) 



id) 



is negative in nematics with positive dielectric anisotropy (Ae = ey — e± > 0) and positive otherwise (/ij| 
0) and fif = //W(6» = tt/2)). This can readily be verified by expanding Eq. (|36|l in powers of the small parameter 



Ae: 



where 



,(d) 



^2 r 



2R\ 



2Ae 



2c, 



1 5(2e s + 1) : 



(e,|+2ej.)/3. 



r P 2 (cos# ) 



(61) 



(62) 



In order to relate anisotropy of to dielectric anisotropy, one needs a relation between the order parameter and 



the dielectric constants. This connection is given by the Maier-Meier theory^ 

3e. 



ey - 1 = 3y 
e± - 1 = 3y 



2e s + l 
3e s 



2c, 



1 



(l + 25 2 F 2 (cos/3)) 
(l-S 2 P 2 (cos/3)). 



(63) 



From Eq. Ae ~ S' 2 P 2 (cos/3). This means that both A/z^? and A/i^ change their sign from negative to positive 
when dielectric anisotropy changes its sign from positive to negative. 
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FIG. 6: The solvation free energy of a dipole solute (in mo/Rf units) vs the angle 6o between the solute dipolar moment and 
the nematic director (Fig. 0; Ro = 1.7a. Solid line refers to en > e± (longitudinal solvent dipole, (3 = 0) and the dashed line 
refers to ey < ej_ (transverse solvent dipole, (3 = 90). Solvent parameters are: rj = 0.47, & = 0.8, y = 0.15, A = 0.3a. At (5 = 0, 
the dielectric anisotropy is positive, ej_ = 1.092 and ey = 2.63; at j3 = 90°, the dielectric anisotropy is negative, e± = 1.746 and 
e || = 1.092. 

The magnitude of solvation anisotropy A/i^ is strongly affected by the size of the solute. The continuum estimate 
of the solvation anisotropy parameter 

K = Avt*/& l& = \(ltf+2$) (64) 

gives very low anisotropies (k ~ 0.02 — 0.03) for a large number of nematics. 2 - Once the dipole correlation effects 

are involved through , anisotropy becomes quite significant, k ~ 0.3S* 2 . This solvation anisotropy results in a 
discontinuity of the solvation chemical potential at the point of the isotropic-nematic phase transition. 
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The significant effect of dipolar correlations on the solvation thermodynamics is seen from comparison of lower 
and upper panels in Fig. [7] which shows steady-state Stokes shift of coumarin-153 dissolved in 4-n-pentyl-4'- 
cyanobiphenyl (5CB) nematogen. This chromophore is widely used as a spectroscopic probe of solvation dynamics 
and thermodynamicsi^SiS Stokes shift dynamics of coumarin-153 in the isotropic phase of a liquid crystalline solvent*^ 
and of coumarin-503SiiS in both isotropic and nematic phases have been reported. The continuum limit in the lower 
panel reveals a much weaker anisotropy in the nematic phase than the full microscopic calculation in the upper panel. 
The continuum calculations are also much lower in the absolute magnitude, which is normally off-set by choosing the 
radius i?o instead of the closest approach distance R\ (Fig. [7] dash-dotted line). However, for the present calculation, 
re-scaling the cavity radius does not fully recover the solvation energy. This result suggests that polar nematics might 
produce stronger solvation than isotropic solvents with a comparable dielectric constant. 

The electrostatic field of coumarin-153 is similar to that of a point dipole^ which makes it a convenient system to 
test our theory. Since the Stokes shift experiments measure only nuclear solvation, the Stokes shift hcAP st (y is the 
wavenumber, cm -1 ) was calculated according to the following expression 

hcAP st = -2//W(m = Am ). (65) 

Here, the difference in the dipole moments in the excited and ground states of the chromophore Atoo is substituted 
for the solute dipole moment. The nuclear component of solvation is calculated in the additive approximation^ in 
which the overall solvation free energy in a polar/polarizable liquid is assumed to be the sum of the nuclear and 
electronic solvation components. The overall solvation free energy /i(e|| , e_L, y c ff) is calculated from the anisotropic 
dielectric constant in the component /i c (e||, e_i_) and the effective dipolar density y c g in fJ-oo(ycs)- The effective dipolar 
density is defined as^ 

y e g = (4TT/9k B T)p(m') 2 + {4w/3)pa, (66) 

where a is the dipolar polarizability and ml is the average dipole moment of the solvent in the liquid. The nuclear 
free energy of solvation is then given by 

Hn = (J,(e\\,e±,yes) - fi{n^,n 2 L ,y e ), (67) 

where n±u is the anisotropic refractive index and the density of induced dipoles is 

y e = (47r/3)ap. (68) 

For coumarin-153, the radius Rq — 4.89 A and the dipole moment change Amo = 7.53 D have been adopted^ 

5CB was chosen as a typical nematogen with its physical properties well documented in the literature: the isotropic- 
nematic transition temperature, Tin = 308.2 K 38 , dipole moment, m — 4.75 Dp22> and the azimuthal angle between the 
dipole moment and the long axis, (3 = 21.6°4£ The temperature dependences of the static dielectric constantS)22ii£ 
refraction indexesj 3 ^^ the order parameter^ and density 3 ^ have also been reported. For the solvent diameter the 
value for benzene, a = 5.27 KM. was adopted, and the polarization correlation length was calculated from Fig. as 
A = 0.35cr corresponding to (to*) 2 = 3.3—3.6. The dipole moment m' was calculated using the Onsager approximation 

, (» 2 + 2)(2e, s + l) 
m = 3(2 C . +n») m ' (69) 

where e s is given by Eq. (|62|l and n 2 = (ra 2 + 2ra 2 L )/3. 

Note that this approximation is made in the Maier-Meier theory [Eq. I|63|l ] which, nevertheless, describes dielectric 
properties of polar nematics reasonably welli^S^ 3 . All results in Fig.0have been obtained at experimentally documented 
parameters of 5CB, the gap between the Stokes shift curves around the isotropic-nematic transition temperature 
reflects the absence of experimental data in this temperature range. 

B. Stokes shift dynamics 

The dynamics of solvation following the photoinduced change in the solute charge distribution is recorded by 
measuring the Stokes shift correlation functionii 
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FIG. 7: Stokes shift of coumarin-153 in nematic and isotropic phases of 5CB. The upper panel shows the microscopic calculation 
according to Eq. 1581 . the lower panel shows the continuum limit. The solid lines refer to the parallel alignment of the solute 
dipole with the nematic director, the dashed lines refer to the perpendicular alignment [Fig. 0. The solute radius and the 
dipole moment are Ro = 4.89 A and Amo = 7.5 D, respectively. The temperature dependent dielectric constants, 22i«i order 
parameter,—- and density— of 5CB are taken from experiment; A = 0.35a. The dash-dotted line in the upper panel refers to 
the continuum isotropic result calculated with the cavity radius equal to Rq. 

where E(t) is the time-dependent energy of the solute. The calculation of this function is normally accomplished 
within the linear response theory. Several formulations of the theory are available in the literature^*^*^ and we 
adopt here the formulation due to Wolynes^ which represents the Laplace transform E(s) of the time-dependent 
function 25(f) as the equilibrium solvation energy characterized by the dielectric constant e(s) 

E(s) = E(e\\(s),ex(s)) = ^ (/x(ey (s), e x (s), y eS ) - n\, j/ e )) • (71) 

The dependence of A on the dielectric constant of a nematic solvent is generally unknown. We will therefore assume 
this parameter independent of the Laplace variable s. Within this approximation, the Stokes shift function is fully 
determined by the continuum expression for the solvation energy [Eq. (|36[) ]. The results for S(t) presented here are 
therefore based on the continuum solution for the dipolar solvation energy in an anisotropic dielectric with axial 
symmetry. 

Experimental dielectric data on nematics are well characterized by the two-exponential form 

£a(s) = £ooa + (eOa - £ooa) V] — , (72) 



where g\ + g2 — 1 and a stands for ||, _L or iso (isotropic). This form reflects two different relaxations: high-frequency 
rotation around a long molecular axis (relaxation time T2) and low- frequency rotation around a short molecular axis 
(relaxation time ti)M> In isotropic phase, T2i SO is about ten times smaller that tu so , and the contribution from high- 
frequency rotation to dielectric loss is usually small and is rarely resolved in the dielectric experiments. 28-38 With 
only one relaxation time tu so , Eq. (|72|l reduces to the Debye dispersion resulting in a single-exponential Stokes shift 
correlation function with the relaxation time ts = [(2eoo + l)/(2eo + 1)]th so fRef. I47T1 . 

The liquid crystalline order of the nematic phase hinders rotations around a short molecular axis when it is per- 
pendicular to the nematic director. The corresponding relaxation time Ty is about ten times larger than tu so . In 
contrary, ru_ is much smaller than tu so , and can be comparable to T2 a (a =||, iso)* — 

With the dielectric constant from Eq. l(72^l substituted into Eq. the function S(t) needs to be calculated 
numerically. Following the procedure described in Ref. lisl the function F(s) = —sE(s) was fitted to a sum of 
Cole-Davidson type functions 

F(s) = (25(0) - iS(oo)) ]T . , (73) 

I 

with the exponents "fi and Stokes shift relaxation times tsv, the linear expansion coefficients are normalized by the 
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condition ^T a,; = 1. Equation (|73|l allows analytical inverse Laplace transform: 

^) = E^^^ r . ( 74 ) 

where r(7, i) is the incomplete gamma function^ 
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FIG. 8: Stokes shift correlation function in the isotropic phase ("iso") and in the nematic phase at two orientations of the 
solute relative to the nematic director. 



Function S(t) calculated for a spherical dipolar dye in 4,4-n-heptyl-cyanopiphenyl (7CB) is presented in Fig. |H1 
Experimental data for dielectric constants, refractive indexes and low-frequency relaxation times is taken from Ref. 
IHol No high-frequency relaxation times T^ a have been reported for this nematogene. S(t) in the isotropic phase (line 
marked "iso" in Fig. |HJ is single-exponential, as expected. 

Two predictions follow from our calculations in the nematic phase. First, S(t) is bi-exponential. The slow component 
with tsi ~ 100 — 500 ns and 71 w 0.82 is related to eii, and the fast component with parameters t$2 ~ 1 — 2 ns and 
72 ~ 1.2 is related to the relaxation of ej_. Second, S(t) is effected by the angle 8q between the solute dipole and the 
nematic director. For nematics with positive dielectric anisotropy (e.g., 7CB) S(t) decays faster for 8 = 90° than 
for 6*o = 0°. Note that the slow component has not been detected in experimental studies of transient Stokes shift in 
nematics (SiiSiii probably because of the limited experimental time resolution (a few ns). All Stokes shift relaxation 
times monotonically increase with lowering temperature in both the nematic and isotropic phase (Fig.[5J). 
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FIG. 9: Temperature dependence of the Stokes shift relaxation times. The upper panel shows the faster relaxation time in 
the isotropic (ts) and nematic (tsz) phases. The lower panel shows the slower relaxation time (tsi) present only below Tin- 
The solid lines refer to the orientation of the solute dipole parallel to the nematic director, while the dashed lines refer to the 
perpendicular orientation. 



17 



V. CONCLUSIONS 



This article presents a microscopic theory of solvation in solvents with axial symmetry. Although applications of 
the theory considered here are limited to nematic liquid crystals, i.e. liquids with inversion symmetry of the polar axis, 
the formalism is also applicable to ferroelectric solvents with a preferential polar direction. This conclusion follows 
from the fact that the linear response approximation leads to a quadratic dependence on the solute dipole moment 
invariant to the dipole flip. For ferroelectrics, however, the solvation energy gains an additional contribution, linear 
in the solute dipole, from the macroscopic polarization of the solvent. 

The full microscopic formulation requires /c-dependent dipolar susceptibility of the nematic solvent which needs 
to be obtained from computer experiment. From the analysis of the results of MC simulations and the microscopic 
formalism, we have derived a formula for the solvation chemical potential which is based on experimentally measurable 
input parameters and a theory parameter, the length of dipolar correlations in the solvent. This correlation length was 
obtained from MC simulations in the range of parameters attainable for the model fluid of dipolar sphero cylinders. 
The understanding of the properties of the correlation length in a broader range of parameters, in particular for higher 
dielectric constants, will require simulations of nematogenes with more realistic intermolecular potentials. 

The theory provides a quantitative framework for interpreting the spectroscopic steady-state and time-resolved 
experiments and makes several experimnetally testable predictions. We show that the equilibrium free energy in the 
nematic phase is a quadratic function of cosine of the angle between the solute dipole and the solvent nematic director. 
The sign of solvation anisotropy is determined by the sign of dielectric anisotropy of the solvent: solvation anisotropy 
is negative in solvents with positive dielectric anisotropy and vice versa. The solvation free energy is discontinuous 
at the point of isotropic-nematic transition. The amplitude of this discontinuity is strongly affected by the size of 
the solute becoming less pronounced for larger solutes. The discontinuity itself and the magnitude of the splitting of 
the solvation chemical potential in the nematic phase are mostly affected by microscopic dipolar correlations in the 
nematic solvent. 

The Stokes shift correlation function in the isotropic phase is one-exponential when dielectric relaxation is given by 
the Debye form. The corresponding relaxation time is a smooth function of temperature through isotropic-nematic 
transition. In the nematic phase, the Stokes shift relaxation becomes bi-exponential with a much slower relaxation 
component related to rotations of polar molecules around their short axes in the nematic potential strongly hindering 
such motions. 

APPENDIX A: DERIVATION OF EQ. (|2?1) 

We start with the formal expression for the solvent dipolar susceptibility^: 
Xs(k) =(m 2 /k B T) / duJidL02 

r - i (A1) 

m(wi)m(w 2 ) S(lji - lu 2 )p(uji) + p(uJi)p(uj2)h(k,LJi,Ld2) , 

where m and lo are the magnitude and orientation of the solvent dipole. Further, the pair correlation function 
/z(k, uj\ , 0J2) is expanded in spherical harmonicsSSiS 

fe(k,u>i,W2) = E E ^(W 1 .ni(w 1 )n.,n a (wa)5 r i:„ 1+ ^(w fc ), (A2) 
where h„„ (k) is the Hankel transform: 



( fc ) = 4 ™ J dr r ii i kr )Klt M • (A3) 


Then the tensor \s can be written as: 

Xs, ni n 2 (k) =^2Xs,n 1 nA k ) Y l*-n 1 -n 2 ( U >k) (A4) 
I 

where coefficients x s ,n 1 n 2 i(k) are proportional to h%f£(k) and depend only on the magnitude of wave-vector k. 
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We next prove the relation 

Uk = 0) - X'(k = 0) = J = 0) (A5) 

where x' is the part of x s which arises from solvent occupied the volume outside of the solute. Since the function 
X' is defined by integrating over the volume twice larger than the solute volume it depends only on the asymptote 
of solvent correlation function on long distances. We will seek the asymptotes of h£ l jL (r) in form 1/r". Since h(12) 
cannot decay slower than the interaction potential, n > 3. According to l|A4Jl and l|A3|l the contribution from the 
asymptotes is proportional to 

nn 

*4W (A6) 



kd 



where d is an arbitrary length larger than size of solvent particle. 
We first consider the case n = 3: 



I = I dx 3 -^ (A7) 



kd 

In the limit k — > 0, this integral converges only if I > 0. Since Xs(k) is analytical function at k = 0, only harmonics 
with I > can have asymptote 1/r 3 . In case of n > 3 we split / into two parts: 

oo 8 

I = fc- 3 / dx^ + fc- 3 / dx&&, (A8) 



fed 

1 1(0.1 -I- IV! 



where parameter 8 is small enough that ji(x) can be replaced by the first term of its Taylor expansion x /(2l + 1) 
Then the first part of / vanishes in limit k — > 0. The second part is proportional to 



s 

l-n+2 



kd 



and in the limit k — + gives non-zero contribution only into term with / = 0. Therefore, in long-wave limit term I = 
m Eq. lHUl arises from harmonics with asymptotes 1/r™ where n > 3, and all therms I > steam from harmonics 
with asymptotes 1/r 3 . Assuming that only the 1/r 3 asymptotes contribute to x'(fc = 0), we arrive in Eq. (| A5|> from 
Eq. I|A4|) and orthogonality of spherical harmonics. 



ACKNOWLEDGMENTS 



This research was supported by the National Science Foundation (CHE-0304694). This is publication #646 from 
the ASU Photosynthesis Center. 



1 M. Born, Z. Phys. 1, 45 (1920). 

2 L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 

3 J. G. Kirkwood, J. Chem. Phys. 2, 351 (1934). 

4 C. J. Cramer and D. G. Truhlar, Chem. Rev. 99, 2161 (1999). 

5 J. Tomasi, Theor. Chem. Acc. 112, 184 (2004). 

6 F. O. Raineri and H. L. Friedman, Adv. Chem. Phys. 107, 81 (1999). 

7 J. Richardi, P. H. Fries, and H. Krinke, J. Chem. Phys. 108, 4079 (1998). 

8 T. Urisu, K. Kajiyama, and Y. Mizushima, Appl. Optics 17, 2366 (1978). 

9 G. Saielli, A. Polimeno, P. L. Nordio, P. Bartolini, M. Ricci, and R. Righini, J. Chem. Soc. Faraday Trans. 94, 121 (1998). 
P. Bartolini, M. Ricci, R. Righini, G. Saielli, A. Polimeno, and P. L. Nordio, Mol. Cryst. Liq. Cryst. 336, 33 (1999). 



19 



11 J. Rau, C. Ferrante, E. Kneuper, F. W. Deeg, and C. Brauchle, J. Phys. Chem. A 105, 5734 (2001). 

12 M. C. B. Mennucci and J. Tomasi, J. Chem. Phys. 102, 6837 (1995). 

13 E. C. B. Mennucci and J. Tomasi, J. Phys. Chem. B 101, 10506 (1997). 

14 M. Inoue and K. Urano, J. Chem. Phys. p. 3672 (1975). 

15 K. Urano and M. Inoue, J. Chem. Phys. 66, 791 (1977). 

16 D. V. Matyushov, J. Chem. Phys. 120, 1375 (2004). 

17 P. de Gennes, The physics of liquid crystals (Clarendon press, Oxford, 1974). 

18 D. Chandler, Phys. Rev. E 48, 2898 (1993). 

19 C. J. F. Bottcher, Theory of Electric Polarization, vol. 1 (Elsevier, 1973). 

20 X. Song, D. Chandler, and R. A. Marcus, J. Phys. Chem. 100, 11954 (1996). 

21 J. P. Hansen and I. R. McDonald, Theory of Simple Liquids (Academic Press, 2003). 

22 C. G. Gray and K. E. Gubbins, Theory of Molecular Liquids (Clarendon Press, Oxford, 1984). 

23 S. H. L. Klapp and G. N. Patey, J. Chem. Phys. 112, 3832 (2000). 

24 J. W. J.M. Caillol and G. Patey, Phys. Rev. A 38, 4772 (1988). 

25 L. Landau and E. Lifshitz, Electrodynamics of continuous media (Pergamon press, Oxford, 1984). 

26 M. Lilichenko and D. V. Matyushov, J. Chem. Phys. 119, 1559 (2003). 

27 S. C. McGrother, A. Gil-Villegas, and G. Jackson, Mol. Phys. 95, 657 (1998). 

28 D. Dunmur, A. Fukuda, and G. Luckhurst, eds., Physical properties of Liquid Crystals: Nematics (Exeter, 2001). 

29 H. E. Stanley, Introduction to phase transitions and critical phenomena (Oxford University Press, New York, 1987). 

30 D. V. Matyushov, Mol. Phys. 79, 795 (1993). 

31 W. Maier and G. Meier, Z. Naturforsch. 16a, 262 (1961). 

32 M. Maroncelli, J. Mol. Liq. 57, 1 (1993). 

33 L. Reynolds, J. A. Gardecki, S. J. V. Frankland, and M. Maroncelli, J. Phys. Chem. 100, 10337 (1996). 

34 P. V. Kumar and M. Maroncelli, J. Chem. Phys. 103, 3038 (1995). 

35 S. Gupta and D. V. Matyushov, J. Phys. Chem. A 108, 2087 (2004). 

36 G. Stell, G. N. Patey, and J. S. H0ye, Adv. Chem. Phys. 18, 183 (1981). 

37 D. V. Matyushov and M. D. Newton, J. Phys. Chem. A 105, 8516 (2001). 

38 S. Urban, B. Gestblom, and A. Wiirflinger, Mol. Cryst. Liq. Cryst. 331, 113 (1999). 

39 S. Urban, B. O. Gestblom, and R. Dabrowski, Phys. Chem. Chem. Phys. 1, 4843 (1999). 

40 J. K. S. Urban and R. Dabrowski, Z. Naturforschg. 55a, 449 (2000). 

41 R. Horn, Le Journal de Physique 39, 105 (1978). 

42 H. Kresse, in Advances in Liquid Crystals, edited by G. H. Brown (Academic Press, 1983), vol. 6, p. 109. 

43 T. K. Bose, B. Campbell, S. Yagihara, and J. Thoen, Phys. Rev. A 36, 5767 (1987). 

44 M. Maroncelli and G. R. Fleming, J. Chem. Phys. 86, 6221 (1987). 

45 P. G. Wolynes, J. Chem. Phys. 86, 5133 (1987). 

46 L. E. Fried and S. Mukamel, J. Chem. Phys. 93, 932 (1990). 

47 B. Bagchi and A. Chandra, Adv. Chem. Phys. 80, 1 (1991). 

48 D. V. Matyushov, J. Chem. Phys. 122, 044502 (2005). 

49 I. S. Gradshteyn and Ryzhik, Table of Integrals, Series, and Products (Academic Press, San Diego, 1994). 

60 M. Davies, R. Moutron, P. A, H, M. S. Beevers, and J. Williams, J. Chem. Soc. Faraday Trans. 2 72, 1447 (1976). 



